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ABSTRACT 

A multi-dimension, time-dependent Monte Carlo code is used to compute sample 7-ray spec- 
tra to explore whether unambiguous constraints could be obtained from 7-ray observations of 
type la supernovae. Both spherical and aspherical geometries are considered and it is shown 
that moderate departures from sphericity can produce viewing-angle effects that are at least 
as significant as those caused by the variation of key parameters in one-dimensional models. 
Thus 7-ray data could in principle carry some geometrical information, and caution should 
be applied when discussing the value of 7-ray data based only on one-dimensional explosion 
models. In light of the limited sensitivity of current 7-ray observatories, the computed theo- 
retical spectra are studied to revisit the issue of whether useful constraints could be obtained 
for moderately nearby objects. The most useful 7-ray measurements are likely to be of the 
light curve and time-dependent hardness ratios, but sensitivity higher than currently available, 
particularly at relatively hard energies (^^ 2 - 3 MeV), is desirable. 
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1 INTRODUCTION 

Although the paradigm that Type la supernovae (SNe la) re- 
sult from the explosions of carbon-oxygen white dwarfs is well 
established, many issues regarding the nature of the progen- 
itors and the explosion mechanism remain unclear (see e.g. 
iHillebrandt & N iemever 2000). Achieving a clearer understanding 
of SNe la is important because of their role in the chemical evolu- 
tion of galaxies and as cosmological distance indicators. 

Considerable theoretical effort has gone into modelling of 
SNe la explosions. Re cently, f ully three-dimensio nal (3D) mod 



elling of the explosion |Reinecke et al...2 002: Game zo et al 
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has become feasible, allowing a detailed treatment of the hydrody- 
namical instabilities and turbulence which play a pivotal role. The 
3D structure predicted by these models has been shown to affect ob- 
servables such a s optical/ultra-violet/infrared (UVOIR) light curve s 
and spectra (e.g. lKasen & Plewal2006l ; ISiml2007l : rSim et alj2007h . 

Although they are primarily detected through their optical 
emission, SNe la are also expected to be 7-ray sources owing 
to the large masses of radioactive i sotopes synthesised in the ex- 
plosion (e.g. iTravaglio et"al] l2004h . In principle, measurements 
of 7-ray emission from SNe la could provide important diag- 
nostics since 7 rays trace almost directly the mass and veloc- 
ity distribution of the products of nuclear burning. Therefore, 
there has be en considerable theoretical work on SN la 7-ray 
emission (e.g.lAm b wani & Sutherland 198i:lBurrows & Thell99oi; 



Miiller et al 



Hoflich et al 



1991; Hoflich etal. 1992; Kuma gai & Nomotdl 19971 : 
19981 ; IGomez-Gomar et al.ll 19981 ; iMilne et al.ll2004h 



Unfortunately, owing to the low sensitivity of 7-ray observatories, 
to date onl y one SN la (SN 1991T) has been detected (see dis- 
cussion by [Milne et^l2004l). However , with current instrumen- 
tation such as the 5 'P/ 7^drenne et alj l2003h on-board Integral 
tWinkler et al .I2OO3I), detection of SN la within several Mpc would 
be possible jOomez-Gomar et al.lll998t) . Moreover, increased sen- 
sitivity in future missions should make detection feasible for more 
distant SNe la. 

Most studies of SN la 7-ray emission focus on making 
predictions from individual models representing specific explo- 
sion mechanisms. In general, the se studies have been r estricted 
to one-dimension (although see iHoflich & Steii] l2002h and to 
Chandrasekhar-mass mo dels (some sub-Chandrasekhar models 
have been considered, e.g. lHoflich et alJl998l ; lG6mez-Gomar et al.l 
Il998l) . Such studies demonstrated that, for nearby SN la, the 
prospects of obtaining useful data are fairly good. Here we adopt 
a different but complementary approach, motivated by the increas- 
ing variety of explosion conditions suggested (both in theoretical 
and semi-empirical studies): rather than determining the quality of 
7-ray data that would be needed to distinguish specific models, we 
attempt to show what might be unambiguously determined solely 
from data and the physics of 7-ray radiation transport. 

We use a multi-dimensional, time-dependent code to com- 
pute 7-ray spectra for a set of parameterised geometries spanning a 
broad range in the relevant physical conditions. Using these refer- 
ence spectra, we highlight the quantities that would be most useful 
diagnostics and therefore most worthy of consideration in the de- 
sign of future instrumentation. 
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In Section|2] we briefly describe the code used to compute the 
reference 7-ray spectra. The reference models we employ are mo- 
tivated in Section [3] and their spectra are discussed in Section |4] 
Guided by our reference spectra, in Section[5]we examine observa- 
tional diagnostics. Finally, in Section[6l we summarise our results. 



2 RADIATIVE TRANSFER CALCULATIONS 

Our calculations were performed using Monte Carlo methods to 
track the emission, propagation, scattering and absorption of quanta 
which represent bundles of 7-ray photons. Monte Carlo methods 
are widely used in modelling 7-ray spectra for SN ejecta, mostly 
owing to the ease with which Compton scattering may be treated in 
a Monte Carlo scheme (see e.g. Milne et al. 2004 for discussion of 
SN la 7-ray transport or iPozdniakov et al..l983. for an extensive re- 
view of Comptonization). They are also highly effective for follow- 
ing radiation transport in SNe la beyond the 7-regime, thereby ob- 
taining bolometric light curve s (following .Cappellaro et al.,,1997.) 
and e ven time-series spectra jMazzali & Lucvlll993 : Kasen et aU 
I2OO6I) . ,— ,— , 

The code used here jSinil2007h is c losely based on the scheme 
and test code developed by Lucvl 1 2005h . We briefly summarise its 
operation and describe the physical processes included. We then 
describe how the histories of the quanta are used to deduce volume 
based emissivities for computation of the emergent 7-ray spectra. 



2.1 Overview of code operation 

As described bv lLucvl j2005h . the Monte Carlo quanta used in our 
code represent indivisible parcels of energy. They begin as pellets 
of radioactive material which were synthesised during nucleosyn- 
thesis in the first few seconds of the SN explosion. 

These pellets decay according to the radioactive nuclei they 
represent and are converted into 7-ray packets. The initial frequen- 
cies of these 7-rays are chosen by sampling the relative probabili- 
ties of the different emission lines that arise from radioactive decay. 

The propagation of the 7-rays is followed in 3D space, time 
and photon frequency until they either escape from the SN or are 
lost from the 7-ray regime; the latter can occur as a result of energy 
loss to non-thermal electrons in Compton scattering or by photoab- 
sorption. During the Monte Carlo simulation, the trajectories of the 
quanta are used to compute volume-based emissivities; the final 
spectra are obtained from a formal solution of the radiative transfer 
equation using these emissivities (see Section[23J. 

Typically, 4 x 10^ quanta are used in each Monte Carlo simu- 
lation. They are tracked on a homologously expanding 100'^ Carte- 
sian grid. The simulations cover the time interval 5 - 150 dy which 
is discretised into 50 logarithmically spaced steps. 

2.2 Radioactive emission 

In the calculations, only ^''Ni and ^®Co are included as sources 
of 7-ray emission. The relev ant nuclear data are taken from 
lAmbwani & Sutherland ( Il988h ; their Table I. All the ^"^Ni is as- 
sumed to be produced in the first seconds of the explosion (i.e. at 
t = 0, for our purposes); ^^Co is produced only as a result of the 
subsequent dec ay of ^''Ni (i.e. no ^^Co i s direct ly sy nthesised). 

Following lAmbwani & Sutherland ( 1 19881) and ILucvI j2005l) . 
positrons released by ^®Co decays are assumed to decelerate 
rapidly and then annihilate in situ, leading to the emission of two 
0.511 MeV 7-ray photons in the local comoving frame (i.e. we 



adopt a positronium fraction of zero; see e.g. lMilneetalj2"004l fora 
discussion). The kinetic energy deposited in the deceleration is as- 
sumed to be radiated at much softer energies and not to contribute 
to the 7-ray spectrum. 

Although other radioactive isotopes are synthesised in the 
early stage of a SN explosion (e.g. **Ti, ^^Ni and ^^Co), they 
are expected to hav e much smaller mass fractions than ^®Ni (see 
e.g. iTravaglio et ^ 12004) and can be neglected at relatively early 
times (i.e. when the light curves are bright and dominated by 
^^Ni and ®^Co). All other sources of 7-ray emission (e.g. nuclear 
excit ation or spa llation) are neglected because, as discussed by 
Gom ez-Gomar et al. ( 1998) , their contributions are negligible until 
very late times. 



2.3 Scattering 

The most important physical process the 7-ray quanta undergo is 
Compton scattering. It is tre a ted vi a the indivisible-energy-packet 
scheme presented by ILucvI ( |2005|) . The Compton cross-section 
((Jc) is given as a function of photon energy and scattering angle 
by the Klein-Nishina formula. In all calculations, we neglect the 
initial electron momentum in the comoving frame. 



2.4 Absorption 

Two 7-ray destruction processes are inclu ded in the ca lculations. 
First, photoabsorption is included foUowing lLucvl ( |2005| ). The pho- 
toabsorption cross-section per unit gram (Kph{E)) is allowed to 
depend on composition in a simple manner: 



Kph{E) = Xpe-grp K28(S) + (1 - Xpe-grp) K.li{E) (I) 

where Xpc-gip is the local mass fraction of iron-group elements 
and K2S (E) and K14 [E) are the photoabsorption cross-sections per 
gram typical of high-and intermediate-mass elements, respectively 
(the subscripts refer to the "typical" atomic mass number, Z). Spe- 
cific e xpressions for k,2s(E) and k,ia{E) are obta ined from Veigel3 
( I1973I) foUowing lAmbwani & Sutherland h988l) . 

Second, for 7-rays above 1.022 MeV, pair production is in- 
cluded. The cross-section for this process, K,pp{E), is allowed to 
depend on composition in the same way as photoabsorption (Equa- 
tion [TJ. The representative cross-sections for intermediate- and 
hig h-mass elements {Z = 14 and 28) are obtained using equation 2 
of jAmbwani & Sutherlan3ll988h . 

When a pair production event occurs, probabilities are as- 
signed to determine whether the outcom e is a positron or a non- 
thermal-electron packet (see ILucvI|2005|) . The kinetic energy of 
both positron and electron are assumed to be lost from the 7- 
ray regime. However, as with the positrons emitted by ^®Co, the 
positrons resulting from pair-creation are assumed to annihilate in 
situ, leading to the emission of two 7-rays at 0.5 1 1 MeV. 



2.5 Extraction of spectra 

As discussed by Lucy (T999', '20031. |2005|) . and in the context of 
the current code by Sim (2007), optimal use of a Monte Carlo 
simulation of radiation transport is generally obtained by record- 
ing volume based estimators and using these to extract the quan- 
tities of interest. The estimators required are those from which we 
obtain the Compton emissivity {r)^ , a function of position, time. 
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frequency and observer line-of-sight direction) and tlie 0.5 11 MeV- 
emissivity arising from pair-production (77*'*'). With these emissivi- 
ties known, a formal solution of the radiative transfer equation can 
be performed to deduce the spectrum (see Lucy 1999; Sim 
The third emissivity term required for the formal solution - describ- 
ing direct 7-ray emission from radioactive decay - does not require 
a Monte Carlo estimator since it can be derived from the assumed 
distribution of radioactive material. 

Monte Carlo estimators for rf" are recorded in every grid cell 
during all relevant time steps in the simulation, on a pre-determined 
frequency grid. In the comoving frame 
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where the summation runs over all 7-ray trajectories inside the grid 
cell during the time step; da /dfi is the Klein-Nishina differential 
cross-section, evaluated at the comoving frequency of the 7-ray 
packet (I'cmf) and for a scattering angle specified by /ic,„f (the co- 
sine of the comoving-frame angle between the direction of packet 
propagation and the observer line-of-sight); v is the ejecta velocity, 
h the rest-frame direction in which the packet is propagating, e,f 
the rest-frame packet energy and ds the trajectory length; / is the 
ratio of photon energies before and after Compton scattering 



/ = 1+^(1- Mcmf) 



(3) 



Ue is the electron number density, V is the volume of the grid cell. 
At the duration of the time step and Ai^ the width of the frequency 
bin. The frequency bin to which a given 7-ray trajectory is assigned 
is determined by Ucmf/f, the comoving frame frequency hypothet- 
ical photons from the 7-packet would have were they to undergo 
Compton scattering into the observer line-of-sight. 

The corresponding estimator for i]^^ is simpler. Since the 
emission is restricted to a single comoving-frame photon energy 
(0.511 MeV), only one estimator per grid cell per time step is re- 
quired for this process; it is given by 
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The summation now runs over all 7-trajectories in the cell during 
the time step and for which the photon energy is above the pair- 
creation threshold. Although 77^,f is a specific emissivity, ry^^j is 
not - hence they differ dimensionally by [Hz^^]. Throughout this 
study it is assumed that emission-line shapes are determined by 
the macroscopic velocity field. Thus the intrinsic emission profile 
associated with rj^^f is taken to be arbitrarily narrow. 

All the spectra shown in this paper were computed via formal 
solutions of the radiative transfer equation using the Monte Carlo 
emissivity estimators given above. The accuracy and validity of the 
use of the estimators was separately tested by comparing spectra 
obtained from the formal solution with those obtained by direct fre- 
quency binning of emergent Monte Carlo quanta for our spherically 
symmetric "control" model (Model SC, see below). To the level of 
the Monte Carlo noise in the frequency-binned spectrum (typically 
a few percent), the two methods produced indistinguishable results, 
as expected. 
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Figure 1. Slices through the x-z velocity-space plane showing the distribu- 
tion of ^"Ni adopted in the models. Grey indicates regions which initially 
contain ^''Ni while black shows the volumes that do not. Panel a shows the 
basic geometry used for Models SC, SS and SFeR. Panels b, c, d, e and f 
show, respectively, Models SM, SFeC, SNiS, AO and AE. All models are 
symmetric under rotation about the z-axis. 



3 MODELS 

We computed 7-ray radiation transport for a set of eight models. 
These are not tied to specific progenitor/explosion scenarios in de- 
tail but are all rooted in previous discussions of SNe la. 

Although quite diverse, our models are not remotely exhaus- 
tive of all possible conditions or geometries. However, as stated in 
Section[T] our intention is not to address specific explosion models 
but to explore whether relatively model-independent information 
can be inferred from the 7-ray regime. 

Compared to other spectral regions (e.g. ultraviolet, optical, 
infrared etc.), the 7-ray spectrum is comparatively simple ~ its 
brightness and shape are mainly determined by only two quanti- 
ties: the instantaneous number of radioactive decays and the opti- 
cal depth to Compton scattering (rc). Thus to first order, 7-ray data 
constrain only the properties that are relevant to these two quanti- 
ties. Of course, other processes (e.g., photoabsorption) have a role 
in establishing the 7-ray spectrum, but their effects are generally 
secondary and they may be disentangled. 

The most easily understood parameter for 7-ray emission is 
the total mass of radioactive material (AfNi). Except for two cases 
(see below), we fix this at 0.6 M0. To zeroth-order, a larger mass 
leads to brighter emission. However, there are additional effects. 

A second relevant quantity is the velocity distribution of ®^Ni 
since Doppler shifts are important in setting emission line widths 
(see e.g. Milne et al. 2004). Detailed radiation transport is not really 
needed to explore this simple relationship, particularly since only 
ext remely sensitive observatio ns could constrain line profiles (see 
e.g. lGomez-Gomar et alll998h . 

Radiation transport is important for the shape of both the 
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light curve and the spectrum. The dominant process affecting 7-ray 
transport in SNe la is Compton scattering. Neglecting light travel 
time effects and relativistic corrections, tc is sensitive only to the 
column density of target electrons (Ne). Since the photon energies 
under consideration (~ 1 MeV) are very large compared to the 
binding energies of electrons in atoms (typically in the 10 eV - 
10 keV range), it is reasonable to regard all electrons as available 
Compton-targets. Thus the Compton optical depth is independent 
of the thermal, ionization and excitation states of the plasma and, 
to a good approximation, also of the density profile along the line- 
of-sight since both the velocity and the velocity gradient appear 
only through Doppler factors. This is in contrast to other wave- 
bands (e.g. the optical) where line opacity dominates giving greater 
significance to the velocity structure. 

In light of these considerations, our models explore a wide 
range of distributions of radioactive source material with tc - as 
a minimal standard of data quality, this distribution should be con- 
strained before data can support or exclude any specific scenario. 

Six of our models are spherically symmetric while the other 
two are aspherical. As mentioned in Section[T] most previous stud- 
ies of 7-ray spectra were restricted to ID calculations. In one of 
the few multi-dimensional studies of SN la 7-ray spectra, Hoflich 
( |2002[) considered an elliptical SN explosion and showed that 
significant viewing angle effects may occur, particularly at rela- 
tively early times. That result is comparable to th ose from stud- 
ies of asp herical SNe la in other wavebands (e.g. Hoflich 199l|; 
iKasen & Pl ewa 20o3; Isim et alj [2007; Hillebrandt et al. 2003) 
which show that viewing angle effects are potentially observable 
and could be confused with other phenomena. Therefore, one of 
our objectives is to add aspherical models to the zoo of other pos- 
sibilities and consider whether they could be distinguished. 

The eight models we use are described below. They are la- 
belled by two to four letters - the first letter identifies whether the 
model is spherical ("S") or aspherical ("A") while the others ab- 
breviate the model properties. The model geometries are shown in 
Fig.[T] All models are assumed to be in homologous expansion and, 
for simplicity, to have uniform mass-density (as noted above, the 
density profile is relatively unimportant for the synthesis of the 7- 
ray spectrum). 

3.1 Model SC: control model 

Our first model, which we use as the standard of comparison for 
all the others, is spherically symmetric with a total mass Mj — 
1.4 M0 and Mm = 0.6 Mq. The maximum ejecta velocity is set 
at Wmax = 1.2 X 10^ cm s~^. This is lower than the maximum 
velocities in real SNe but representative of the velocity regime in 
which the density is large and is therefore appropriate given our 
choice of a uniform density profile (see above). In this model, the 
®^Ni is located in a spherical region at the centre of the ejecta. The 
volume of this region is fixed by an adopted ^^Ni mass fraction of 
0.9. The region outside the central concentration of Ni is assumed 
to be composed of intermediate mass elements. 

3.2 Model SS: super-Chandrasekhar mass model 

Model SS has the same geometry and Ni-distribution as Model SC 
but its densities are everywhere higher by a factor of 1.5. Thus, 
this model is "super-Chandrasekhar", having Mt = 2.1 Mq and 
Mm = 0.9 Mq. This model explores 7-ray emission from ob- 
jects in which tc is higher than in Model SC. Although out- 



side the standard SN la paradigm, Super-Chandrasekhar explo- 
sions have been proposed f or unusually bright events suc h as 
SN 199 IT (Fisher etal."l999'), SN 2003fg dHowell et alj2006h and 
SN 2006gz fflicken et al.ii2007.) . 

3.3 Model SFeR: Fe-rich ejecta 

This model isolates the effect of composition. It adopts the same 
geometry, mass and Ni-distribution as Model SC. However, it as- 
sumes that the majority of the material that is not initially ^^Ni con- 
sists of other heavy nuclei; the photoabsorption and pair-production 
cross-sections are thus different from Model SC. 

3.4 Model SM: well-mixed model 

This model is identical to Model SC (Mt = 1.4 Mq, Mm = 
0.6 Mq) except that ^^Ni is distributed uniformly throughout the 
volume of the ejecta (see panel b of Fig.[I). 

3.5 Model SFeC: stable Fe core 

Recentlv.^ Mazzali et al.l ( 1200 7^ have argued for a model which may 
account for most SNe la. Analysing spectra of a large sample of 
SNe, they concluded that the characteristic structure consists of a 
low-velocity core of stable iron-group material surrounded by a 
^^Ni-rich region. Outside this region, the material underwent in- 
complete nuclear burning and is dominated by intermediate-mass 
nuclei. Motivated by their study, we include a model based on such 
a structure (the geometry is shown in panel c of Fig.[Tll. The inner 
stable iron-group core has a mass of 0.2 Mq. A pure ^*'Ni-region 
lies immediately above this core and contains 0.6 Mq. The outer- 
most region is composed of intermediate-mass elements only. 

3.6 Model SNiS: Ni-rich core and surface layer 

This model differs from all others as it includes two physi- 
cally separated ^"Ni-rich regions: a massive core and a thin sur- 
face layer. This geometry provides an alternative to the "super- 
Chandrasekhar" scenario (see Section 3.2) for unusually bright 
events. For this model we retain Mt — 1.4 Mq, adopt a core 
^''Ni-mass of 0.9 Mq and assume that the outermost 0.1 Mq of the 
ejecta is also pure ^®Ni. This outer layer is motivated by the early - 
time spectroscopic behaviour of SN1991T jMazzali etalj|l995h . 
The material sandwiched between the core and the outer Ni-rich 
shell is taken to be composed of intermediate mass elements. This 
model has the largest Aim in our study, as required to account for 
the brightness of events such as SN1991T. 

3.7 Model AO: off-set Ni 

Here, the basic model is identical to Model SC but the centre-of- 
mass of the Ni ball is displaced along the z-axis by 10 per cent 
of the outer ejecta radius. T his leads to an as pherical model, very 
similar to those discussed bv lSimetalJ ( l2007h . 

3.8 Model AE: ellipsoidal ejecta 

In this model, the maximum ejecta velocity varies with direc- 
tion, producing an el lipsoidal SN. T his is the same basic geome- 
try as considered by iHoflicri ( l2002h . Ellipsoidal explosions have 
been suggested following detection of polarisation in SNe la (e.g. 
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Figure 2. Distribution of source """Ni witli electron scattering optical depth 
(computed in the Thompson Umit) for the models. The distribution is nor- 
malised to drc = 1. The symbols are indicated in the figure. 
Model SFeR is not plotted since, with respect to this distribution, it is iden- 
tical to Model SC. For both the aspherical models (Models AO and AE), 
two curves are plotted - these represent the extremes with respect to view- 
ing angle: for Model AO, the solid line is for a line-of-sight parallel to the 
direction in which the Ni-blob is displaced from the centre of the SN while 
the dashed line is for the anti-parallel light-of-sight; for Model AE, the solid 
line is for viewing down the short axis of the ellipsoid while the dashed line 
is for a line-of-sight along the long axis. 



iHowell et all 1200 ll : Imng et alj|2003h . The maximum velocity ( 
1.5 X 10^ cm s"'^) occurs at the equator. The axis ratio is set to 
5:4; therefore the terminal velocity of the ejecta in the polar di- 
rection is 1.2 X 10^ cm s~^. The Ni-distribution is ellipsoidal in 
the same sense as the total mass distribution and remains centrally 
concentrated. As in all the other models, the mass-density remains 
uniform. 
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3.9 Distribution of emission with opacity 

As discussed above, since the 7-ray spectrum is primarily deter- 
mined by Compton scattering, we expect that a highly relevant 
quantity is the distribution of ^''Ni with re (see Fig.O. tc is a de- 
creasing function of both time and photon energy. However, since 
the energy dependence is universal and, in an homologous flow, 
the time dependence is very simple (tc oc t~^), it is sufficient to 
compare the distribution for one photon energy at one time. For 
convenience, the optical depths used in Fig. [2] are computed using 
the Thompson limit to ac at time t = 50 dy. 

Fig.|2]shows that the models considered are quite diverse and 
reasonably complete, covering a range of plausible single-peaked 
distributions. At the epoch considered, the radioactive material in 
Model SC is concentrated at moderate optical depths, and its dis- 
tribution is bracketed by the extremes shown in Model AO. Mod- 
els SFeC and AE have lower optical depths, but their distributions 
are still peaked significantly away from rc = 0. Model SM has an 
extreme distribution, with a significant fraction of the source ma- 
terial at very low optical depths (rc ~ 0). At the opposite limit. 
Model SS has higher opacities, resulting from the larger density. 
Model SNiS differs from the others in having two peaks in its dis- 
tribution - most of the Ni lies at moderate optical depths but the thin 
surface layer provides a distinct peak at very low optical depth. 

Model SFeR is not shown in Fig. [2j its distribution of ''''Ni 



Figure 3. Representative time series of 7-ray spectra computed with four 
spherically symmetric models (identified in the second panel) at times after 
explosion as indicated on each panel. The ^''Ni and ^^Co emission lines are 
identified in the top panel. The four flux-bands discussed in Section 5.3 are 
indicated in the bottom panel. Fluxes are for a source distance of 1 Mpc. 

with TC is identical to Model SC so that it can isolate the signatures 
of photoabsorption. 



4 DISCUSSION OF SPECTRA 

Fig.|3]shows a time series of spectra computed from Model SC and 
spectra from three other spherical models (SS, SFeR and SM) are 
over-plotted. (For clarity, spectra from the remaining four models 
are not plotted. They are, however, included in all the discussions 
of observable diagnostics in subsequent sections.) 

At times close to maximum light the 7-ray spectrum consists 
of strong emission lines, mainly due to ''''Co, with significant con- 
tinuum arising from Compton scattering of line photons. Since op- 
tical depths decrease with time in the expanding ejecta, the strength 
of the lines relative to the continuum increases with time. At early 
times, the "mixed" model (Model SM) shows a harder continuum 
than the standard model and stronger Ni emission lines (particularly 
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0.158 MeV and 0.275 MeV). These lines are also fairly strong in 
Model SNiS (not pl otted) where they originate in the surface layer 
of ^'^Ni. As noted bv lGomez-Gomar et al.] l fl99i) . these Ni lines can 
only form when the source Ni lies at small rc ; since oc decreases 
with increasing photon energy, the soft-energy lines are degraded 
most easily, becoming swamped by photons down-scattered from 
harder energies. 

In Model SS, the 0.158 and 0.275 MeV lines are almost com- 
pletely buried by the strong continuum which persists until well af- 
ter maximum light as a consequence of the high re's in this model. 

Model SFeR shows the effect of composition as intended. 
Above about 0.3 MeV, its spectra are indistinguishable from those 
of Model SC. At soft energies, however, the Model SFeR flux is 
lower by up to an order of magnitude. This is due entirely to the 
difference in the photoabsorption cross-section arising from the 
choice of composition. Although weaker, a similar effect appears in 
Model SNiS because of the relatively large amounts of iron-group 
material in that model, particularly in the surface layer. 

Should energy-resolved data of sufficient sensitivity to mea- 
sure both the line and continuum emission across the entire 0. 1 - 
3.5 MeV spectrum be obtained, it would be used to evaluate SN 
models by direct comparison. However, given the rarity of very 
nearby SNe la and the limited sensitivity of 7-ray observatories, 
in the next section we will consider what information could be ex- 
tracted from low quality 7-ray data in the form of simple diagnos- 
tics. We will focus on line and continuum fluxes. Although the en- 
ergy resolution of 7-ray instruments ca n be high enough to resolve 
spectral lines (e.g. with SPI Integral, see lRogues et al.l2003h . sensi- 
tivity limits make it practically impossible to measure detailed line 
shape s except for extremely nearby events dOomez-Gomar et"al] 
Il998h . 



5 OBSERVATIONAL DIAGNOSTICS 

We consider three different potential diagnostics that could, in 
princple, be obtained from 7-ray data. Throughout our discussion, 
we use theoretical source count rates to illustrate the value of the di- 
agnostics for distinguishing our models. We do not, however, con- 
sider any additional effects (e.g. instrumental noise) which would 
make detections much har der with current observatories (see e.g. 
iGomez-Gomar et alj|l998h . 

We begin with the simplest diagnostic, the energy-integrated 
light curve. 



5.1 7-ray light curves 

5.1.1 Discussion 

The 7-ray light curve is the most useful for constraining Mm; since 
in this spectral region the decay products of ^®Ni and its daughter 
nucleus ^®Co are seen most directly. 

Fig. m shows energy-integrated (0.1 - 3.5 MeV) light curves 
from the models. All the light curves have a characteristic, single- 
peaked shape. They peak later (tpcak = 50-80 dy) and are simpler 
in shape than model bolometric (or optical/infrared) light curves 
(see e.g. B linnikov et al. 2006; Kasen et al. 2006; Kasen 2006; 
iKasen & Wooslevil2007l) since photon trapping plays a far less sig- 
nificant role in the 7-ray regime. 

After about 100 dy, all our models are fairly optically thin and 
the 7-ray luminosity is an unambiguous measure of A/ni via 
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Figure 4. Model light curves for the 0.1 - 3.5 MeV 7-ray region. The 
flux is given for a source distance of 1 Mpc. Times are measured from the 
explosion. For the aspherical models (AO and AE), two curves are plotted 
- as in Fig. 12] these represent the extremes vi'ith respect to viewing angle. 



L-f{t> 100 dy) X ecoMNiexp(-f/tco)Aco 

^ 1.23 X 10''=*^ exp(-t/tco) ergs s^^ 

(5) 

where L~, (t) is the emitted 7-ray luminosity, eco is the total 7-ray 
energy released during the decay of one gram of ■"^''Co and tcc IS 
the ^''Co lifetime. In Fig. [4] this optically-thin behaviour causes the 
light curves of all models with the same adopted AImi (SC, SFeR, 
SM, SFeC, AO and AE) to converge for t > 100 dy while models 
with higher Mni (SS and SNiS) are proportionally brighter. 

At earlier times, individual model light curves differ signif- 
icantly - the spread is about a factor of two at t ~ 40 dy. The 
light curves of models containing Ni at relatively low rc (SM, AE, 
SNiS, SFeC and face-on AO; see Fig.[2ll tend to rise quickly and 
peak early but the light curve shapes are different and do not ad- 
here to a one-parameter family. The aspherical models show view- 
ing angle dependence of up to 50 per cent at early times, decreasing 
to about 10-20 per cent at tpcak which is comparable to differences 
between our various spherically symmetric models. 

The integrated light curve is almost insensitivity to photoab- 
sorption (in Fig. |4] compare Models SFeR and SC) - most of the 
7-rays are too hard for significant photoabsorption to occur, regard- 
less of the composition of the absorbing material. 

5.1.2 Diagnostic value 

With good time coverage (40 - 120 dy after explosion), the 7-ray 
light curve alone can provide useful information: at late times it 
gives the most direct measurement of Mm while at earlier stages 
it is sensitive to the distribution of rc (i.e. electrons and therefore 
total mass). 

To assess the diagnostic value of a SN la 7-ray light curve, we 
converted our theoretical light curves to source counts. For illustra- 
tive purpos es, we used the energy-dependent effective area of In - 
tegral SPI ( Atti e et alj2003l;lDietil et alj2003l ; [srumer et al.l2003l) . 
adopted a distance of 18 Mpc to the SN e (an approximate distance 
to the Virgo cluster. lFougue et al.l200ll) and assumed an integration 
time of 10^ s (~ 1 dy). Under these conditions Model SC would 
give ~ 460 photon counts at t ~ 50 dy, decreasing to ~ 310 counts 
at t ~ 100 dy. 
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Therefore, if photon counting statistics were the only hmita- 
tion, it would be quite feasible to deduce SN la ''''Ni-masses to a 
precision of about 10 per cent for objects at distances up to Virgo. 
Some rough information on the light curve shape between about 40 
and 100 dy after the explosion could also be obtained. The level of 
precision would be comparable to the typical differences between 
the models investigated here, so that a time sequence of data could 
place fairly robust constraints not only on the total ^''Ni mass, but 
also on its distribution with opacity. 

Of course, these estimates are very crude and optimistic since 
additional sources of error (e.g. background) have not been consid- 
ered, but they do suggest that, in principle, useful information is 
contained in light curve data. 

5.2 Line ratios 

5.2.1 Discussion 

7-ray emission line peak-inte n sity r atios provide sim- 
ple diagnostics teurrows & Thd Il990l : iHoflich et alj 1 19981 : 
iGomez-Gomar et alj|l998l) . As in earlier studies, we examine only 
a few lines and their ratios in detail. Owing to the simplicity of 
spectral formation in the 7-ray regime, exactly which lines are 
considered is fairly unimportant - the qualitative behaviour of any 
line ratio will be similar to that of one of those we show. 

Since the strongest lines originate in the decay of ^®Co, their 
intrinsic relative strengths are set by nuclear physics alone - thus 
late-time measurements of line ratios yield no information that can- 
not be gleaned from the energy-integrated light curve. However, at 
early times, the energy sensitivity of ac ma kes the intensity ratio 
of lines at different energies sensitive to A^e ^Hoflich etal.lll998h . 

The left panel of Fig. [5] shows the ratio = F(0.847 
MeV) /i^(2.598 MeV) . Both lines in this ratio are from the decay of 
^®Co. We chose to consider this particular pair of lines since they 
are strong, unblended and well-separated in energy. Although an 
even more widely spaced pair of lines associated with ''''Co could 
be achieved considering the 0.5 1 1 keV positron-annihilation line, 
we prefer to avoid reliance on our adopted positronium fraction. 
The R\ ratio in Fig.[5]has been normalised such that it tends asymp- 
totically to 1 .0 in the optically thin limit. 

For times up to about 70 dy, R\ discriminates between dif- 
ferent distributions of ^"Ni with tc (see Fig. |2](. Using this ratio, 
spherical models in which the radioactive material is behind most 
opacity (Model SS) and least opacity (Model SM) can be easily 
discriminated from the control model. Aspherical models behave 
similarly - in Model AO, Ri is largest if observed from the side to 
which the Ni blob is displaced; the effect is present but weaker in 
Model AE because tc is typically smaller. 

As for Ri , for any other pair of strong Co emission lines the 
flux ratio F{e\) / F[e2) will be determined by the ratio of Comp- 
ton opacity at the two line energies, crc(e2)/o"c(ei). We will not 
discuss these other potential diagnostic ratios but only comment 
that the most useful ratios will always be those between lines of 
significantly different energy (as in 7?i). Thus, the ratio of the two 
strongest lines (F(0.847 MeV)/f'(1.238 MeV)), although easier 
to measure, has less diagnostic value than R\ . 

The second line-flux ratio shown in Fig.|5]is between the ^"Ni 
0.158 MeV line and the ^''Co 0.847 MeV line, referred to as R2. In 
the figure, it is normalised to the optically thin limit, a decreasing 
function of time set by the decay rates of ^®Ni and ^''Co. 

For most models, R2 gives the same relative ordering as R\ at 
early times. Again, this is primarily due to the energy-dependence 



of ac. However, the effect is weaker in R2 than R\ because it 
is opposed by an increase in the Compton continuum flux around 
0.158 MeV. Disentangling these two effects would require suffi- 
ciently high quality data that the 7-ray continuum level around the 
0.158 MeV line can be measured. 

Unlike R\, however, R2 effectively separates Model SFeR 
from any of the others at t < 50 dy since the 0.158 MeV line is 
sufficientl y soft to be affected b y photoabsorption. Thus, in agree- 
ment with lGomez-Gomar et alj ( 11998) , we conclude that R2 (or an 
equivalent ratio of one of the low-energy Ni lines to a harder Co 
line) provides an important diagnostic for composition. 



5.2.2 Diagnostic value 

The Ri and R2 ratios defined above provide useful diagnostics for 
two of the quantities that may be constrained via the 7-ray spec- 
trum, the distribution of total mass and the composition of the pho- 
toabsorbing plasma. 

Peak count rates for the relevant emission lines computed 
adopting the same conditions described in Section [5.1.2| are shown 
in Fig.|6] For both the 0.158 and 0.847 MeV lines, the peak count 
rates are several thousand/Me V. The intrinsic line widths, deter- 
mined by the velocity range in the ejecta, are ~ 6 keV for the 
0.158 MeV line and ~ 34 keV for the 0.847 MeV line. These 
widths are large enough to be moderately well-resolved by mod- 
em instruments ( Attie et al. 2003; Roques et al. 2003). For our as- 
sumed distance and integration time, one would only expect around 
50 source counts in the 0.847 MeV line and fewer than 10 source 
counts in the 0.158 MeV line. 

Owing to the combination of fewer source photons and 
smaller adopted detector effective area, there are fewer counts/MeV 
for the 2.598 MeV line - only a handful of counts integrated over 
the line profile. 

Thus, statistical errors caused by small number counts are very 
substantial in both R\ and R2, meaning that these ratios could not 
contain useful information for objects at our chosen distance. A 
tenfold increase in source counts would be needed to make the sta- 
tistical accuracy comparable to any of the differences between the 
models (as shown in Fig.|5}. With 50 - 100 times more counts and 
good time coverage these ratios could contain detailed information 
although additional instrumental effects would further limit their 
value in practise. 

In light of this difficulty, in the next section we consider an 
alternative approach in which the spectrum is coarsely binned and 
hardness ratios are extracted. 



5.3 Hardness ratios 

5.3.1 Discussion 

Given the simplicity of the 7-ray spectrum, moderately complete 
information can be obtained from hardness ratios alone; this is in 
contrast to other wavebands where spectra are complex and not 
well-described by photometry alone. Some prospects of using rel- 
atively broad energy bands for pr obing the soft, continuum re gions 
of the spectrum were discussed bv lGomez-Gomar et (11998). We 
extend this to consider hardness ratios involving both continuum- 
dominated and higher energy line-dominated bands. 

We have divided the spectrum into four energy bands, two in 
which the Compton continuum contributes significantly to the flux 
(Ci and C2) and two in which strong line emission always domi- 
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Figure 5. Peak line flux ratios from the various models as functions of time since explosion. The left panel shows the ratio R\ = P'(0.847 MeV)/i<'(2.598 
MeV) and the right panel R2 = F(0.158 MeV)/F(0.847 MeV). For the two aspherical models (Models AO and AE), two curves are plotted - as in Fig.|2] 
these represent the extremes with respect to viewing angle. Flux ratios are normalised to the optically thin limit (thus, the normalised R\ -ratio asymptotes to 
1.0 while R2 asymptotically approaches the instantaneous mass ratio of ^^Co to ^^Ni). 




Figure 6. Predicted peak counts MeV"! for the 0.158, 0.847 and 2.598 MeV emission lines computed adopting the Integral SPI effective area, a recording 
interval of 10^ s and a source distance of 18 Mpc. For both the aspherical models (Models AO and AE), two curves are plotted - as in Fig. [2] these represent 
the extremes with respect to viewing angle. 




Figure 7. Hardness ratios computed from the various models as functions of time since explosion. For the aspherical models, two curves are plotted - as in 
Fig. [2] these represent the results for the two extreme observer lines-of-sight. 
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Figure 8. Synthetic ligiit curves for the various models in each of the four energy bands considered computed adopting the Integral SPI effective area, a 
recording interval of 10^ s and a source distance of 18 Mpc. For the aspherical models, two curves are plotted - as in Fig.|2] these represent the results for the 
two extreme observer lines-of-sight. 



nates (Li and L2). The bands are defined below and are indicated 
in Fig. [3] 

• Ci (0.1 - 0.3 MeV): emission is generally dominated by 
Compton down-scattering of harder photons, photoabsorption can 
be significant. 

• C2 (0.3 - 0.7 MeV): emission is predominantly through 
Compton down-scattering. There is some direct line emission, 
mostly via the 0.51 1 MeV electron-positron annihilation line. 

• Li (0.7 - 2 MeV): emission is dominated by strong lines of 
Co supplemented by a moderately weak Compton continuum. 

• 1/2 (> 2 MeV): lines of Co again dominate, now with an even 
weaker Compton continuum. 

The line ratios Ri and R2 (see Section [5!2l ), are affected by 
the change in ac between different energies. A similar effect can 
be anticipated in the hardness ratio between the two line-dominated 
bands (Li and I/2) - such an effect is present but weak (see below). 
Another effect, however, can be noted when comparing fluxes at 
soft energies. The strength of the Compton continuum relative to 
the lines increases with opacity and thus high optical depths lead to 
softer hardness ratios when Ci or C2 are involved. 

This is illustrated in Fig. [7] which shows four hardness ratios 
computed with different combinations of the four 7-ray bands. The 
hardness ratios H\ = Ci/Li, H2 = C2/L1 and H3 = C1/C2, all 
show the effect of high tc in softening the spectrum. Of these. Hi 
and H2 are best at distinguishing the models, although the scale of 
the differences is smaller than in the line ratios (c.f. Fig.|5l(. 

The hardness ratios involving the Ci band (Hi and H3) are 
sensitive to the composition of the scattering material - these ratio 
clearly separate Model SFeR from Model SC. This behaviour is 
analogous to that of the R2 line ratio. 

The fourth ratio, H4 = L1/L2 shows an opposite trend - it is 
hardest when optical depths are high (Fig.[5](. As mentioned above, 
this can be understood by analogy to Ri but H4, is less sensitive 
since it represents an average over a range of energies. 



Thus, 7-ray hardness ratios can contain quantitative informa- 
tion similar to that available from the line ratios discussed in Sec- 
tion l5.2l Although hardness ratios are less sensitive than line ratios, 
they may benefit from increased accuracy owing to the broad en- 
ergy ranges over which the spectra can be integrated. 

5.3.2 Diagnostic Value 

Light curves of the various models in each of the four energy bands 
defined in Section 15.3.1 1 were computed, again adopting the same 
conditions described in Section [5.1.2l and are shown in Fig.[8l 

In the three softer bands, we expect roughly 100 - 150 source 
counts from most models at ~ 40 dy with the adopted observation 
parameters, giving an optimistic estimate of ~ 15 per cent statisti- 
cal error in Hi, H2 and Hg. This is sufficient to separate the most 
different models - Models SC and SFeR could be distinguished 
via either Hi or H3 while Models SM and SS would be measur- 
ably different from each other and marginally distinguishable from 
Model SC. However, there is little prospect of using these ratios to 
find any of the effects of asphericity captured by Models AO and 
AE. 

The predicted L2-band source count rate is too low to allow 
the H^-mtio to be useful - only a handful of counts are anticipated 
in 1/2 for any of the models, a consequence of the reduction in both 
our adopted detector effective area and number of source photons 
at hard energies. 

Still, hardness ratios could carry some diagnostic information. 
In contrast to the line-ratio diagnostics (Section [5.2.2| l, an increase 
in the count number by even a modest factor would allow fairly 
detailed analysis. An increase by a factor of ~ 10 - say as a result of 
a serendipitously nearby SN la - would reduce the source statistical 
error in the Hi, H2 and H3 ratios to less than 10 per cent. At this 
level of accuracy, these ratios would contain quantitatively useful 
information, particularly if their time evolution could be followed 
through the peak of 7-ray emission. 
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6 SUMMARY AND PROSPECTS 

Using a Monte Carlo code we computed 7-ray spectra for a vari- 
ety of models to explore whether unambiguous constraints could 
be obtained from 7-ray observations of SNe la. Two aspherical 
toy geometries (a lop-sided distribution of Ni and an ellipsoidal 
ejecta) show that moderate departures from sphericity can produce 
viewing-angle effects at least as significant as those due to varia- 
tions of key parameters in ID models. Thus 7-ray data could carry 
some useful constraints on possible geometries, but caution must 
be applied when evaluating its potential usefulness in distinguish- 
ing specific explosion scenarios. 

Given the limited sensitivity of current 7-ray mis- 
sions, we conclude , in a greement with previous studies (e.g. 
lOomez-Gomar et al] Il998 "). that there is little prospect for ob- 
taining useful constraints from line-ratio diagnostics except for 
fortuitously nearby objects. Instead, we suggest that the best 
prospects are offered by hardness ratios. In particular, owing to the 
simplicity of the physics underlying the 7-ray spectrum, a simple 
ratio of the total emission in a hard energy, line-dominated part 
of the spectrum to a soft energy, continuum-dominated region 
could discriminate between our more extreme models. We also 
emphasise the value of obtaining multiple observations over a 
wide time period given the diagnostic power of the light curve 
shape. In planning future 7-ray missions, greater sensitivity at 
harder energies ( ^ 2 MeV) should be given high priority since 
several of the most potentially useful diagnostics (e.g. our Ri - and 
_ff4 -ratios) require measuring the hardest energies in the spectrum. 
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